The cos 20 azimuthal asymmetry of unpolarized pp collisions at Tevatron 
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We calculate the cos 2(j> azimuthal asymmetry of the unpolarized pp Drell-Yan dilepton production 
process in the Z-resonance region at the Tevatron kinematic domain. Such an azimuthal asymmetry 
can provide additional information about a spin-related new parton distribution function, i.e., the 
Boer-Mulders function of the proton, compared to the pp process. Therefore the available data of 
unpolarized proton-antiproton collision at Tevatron can contribute to our study on the spin structure 
of the nucleon. 
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The study of the intrinsic transverse momentum de- 
pendent (TMD) distribution functions has received much 
attention in recent years [ij. Such new quantities 
of the nucleon provide us a significant perspective on 
understanding the spin structure of hadrons and the 
non-perturbative properties of quantum chromodynam- 
ics (QCD). The intrinsic transverse momentum of par- 
tons may cause special effects in high energy scatter- 
ing experiments [3]. It was naively speculated that the 
polarization of at least one incoming hadron is neces- 
sary to investigate the spin-related structure and proper- 
ties of hadrons, however the situation will change if the 
transversal motions of quarks inside the hadron will take 
into account. The Drell-Yan process is an ideal ground 
for testing perturbative QCD and probing TMD distri- 
bution functions, and its cross section is well described by 
next-to- leading order QCD calculations [3] . Surprisingly, 
the first measurement of the Drell-Yan angular distribu- 
tion, performed by NAlO Collaboration for ttN at 140, 
194 and 286 GeV, indicates a sizable cos2(/) azimuthal 
asymmetry 0, 0| which cannot be described byleading 
and next-to-leading order perturbative QCD [u\. Fur- 
thermore, the subsequent result by the Fermilab E615 
Collaboration reveals that the Lam- Tung relation Q, 
which is analogous to the Callan-Gross relation in 
deep-inelastic scattering, obtained as a consequence of 
the spin-i nature of the quarks, is clearly violated 
The violation has also been tested in recent pd and pp 
Drell-Yan dimuon processes measured by E866/NuSea 
Collaboration l^, 11|. 



Several attempts were made to interpret this asym- 
metry, such as the factorization breaking QCD vacuum 
effect [Qlwhich is possible the helicity flip in the instanton 
model [13, higher twist effect |13l - ll5| and the coherent 
states [16[. Boer pointed out that the cos 20 azimuthal 
asymmetry could be due to a non-vanished TMD dis- 
tribution function hi{x,p'^) [l3|, named as the Boer- 
Mulders function later, as one of the eight leading-twist 



TMD distribution function contained in [18l US 
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where $ is the quark-quark correlation matrix, defined 
as 
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The Boer-Mulders function is another time-reversal odd 
(T-odd) distribution function which characterizes the 
correlation between quark transverse momentum and 
quark transverse spin, analogous to the Sivers function 
f iri^ J Pt) which signifies the correlation between quark 
transverse momentum and hadron transverse spin [20| . 
The non-vanished T-odd distribution functions can arise 
from the initial-state or final-state interaction ^21r.-24j. 
In general, the path-order Wilson line arising from the 
requirement of QCD gauge invariance for quark correla- 
tion functions provides non-trivial phases and le ads to 
non- vanished T-odd distribution functions 
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to the present of Wilson line, opposite sign of the Boer- 
Mulders function or Sivers function in semi-inclusive deep 
inelastic scattering (SIDIS) and Drell-Yan processes is 
expected 2^ 30| . 
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The existence of T-odd distribution function can 
cause azimuthal asymmetries in SIDIS at leading twist 
level [l^, and the product of two Boer-Mulders func- 
tions of two incoming hadrons may give a sizable cos 20 
azimuthal asymmetry in unpolarized Drell-Yan processes 
by establishing a preferred transverse momentum di- 
rection from the spin-transverse momentum correlation. 
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which is called the Boer-Mulders effect [17|. Thus, the 
measurement of the Boer-Mulders function will promote 
our understanding of QCD. Many theoretical and phe- 
nomenological studies are carried out along this direc- 
tion IH-El. 

Recently, the Collider Detector at Fermilab (CDF) 
Collaboration first measured the angular distribution co- 
efficients of Drell-Yan e+e~ pairs in the Z mass region 
from unpolarized pp collisions p -t- p — > 7*/Z + X — > 
Z+r-HXatys=1.96 TeV [ll]. This indicates that 
it is feasible to investigate spin physics at Tevatron. In 
this paper, we calculate the cos 2(/) azimuthal asymmetry 
caused by the Boer-Mulders effect in the Z-pole region 
with the kinematic conditions at Tevatron. 

The angular distribution coefficients are generally 
frame dependent. We choose the Collins-Soper (CS) 
frame [50| to perform the calculation. It is the center 
of mass of the lepton pair with the z axis defined as the 
bisector of p and p beams. The polar angular Q is defined 
as the angular of the positive lepton with respect to the 
z axis direction, and the azimuthal angular is defined 
as the angular of the lepton plane with respect to the 
proton plane. In this frame the Lam- Tung relation is in- 
sensitive to the higher fixed-order perturbative QCD [Flj 
or the QCD resummation 52l454|. The angular differen- 



tial cross section for unpolarized Drell-Yan process has 
the general form: 



the fine structure constant, the mass of proton, the trans- 
verse momentum, and invariant mass of ^* /Z respec- 
tively. The structure function notation in this equation 
is defined as 
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where px, kx are the transverse momenta of quarks in 
proton and antiproton, and h = ^ is the direction of 
the transverse momentum of ^* /Z . The coefficients Ki, 
K3 and K4 are expressed as: 
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where Ca is the charge of quarks (antiquarks) , and gv and 
gA are the vector and axial-vector couplingconstants to 
the Z-boson. We take their values in Ref. [56!]. The Ci is 
defined as: 
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where $7 is the solid angle and A, /i, and v are angu- 
lar distribution coefficients. For azimuthal symmetrical 
scattering, the coefficients 11 = v = {). It can also be 
written as 0, HI] : 

^ =Wt(1 + cos^ 61) -t- Wl{1 - cos^ 6) 
dil 

+ Wa sin 29 cos </) -I- H/aa sin^ 9 cos 20. (5) 

When taking into account both virtual photon and Z- 
boson contribution, the leading order unpolarized Drell- 
Yan cross section is \vA 
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where xi, X2 are the Bjorken variables standing for the 
longitudinal momentum fractions carried by the partons 
in the proton and antiproton, and a, M, qt, and Q are 



where j = ^ or a. The Z-boson propagator Xi is given 
by: 
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where 9w is the Weinberg angle. In Eq.®, we as- 
sume that the TMD distribution functions for anti- 
quarks (quarks) in the antiproton are the same as those 
for quarks (antiquarks) in proton, and the summation 
over the index a are for different flavors with a = u, d, u, 
and d. 

In our calculation, we take the Boer-Mulders functions 
extracted from pD and pp Drell-Yan data 43]. The 
parametrizations for hi^{x) is (isj : 

hi^{x)^H,x-^{l-x)''h,{x), (15) 
and the TMD part is parametrized with a Gaussian form: 
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This parametrization is based on the assumption that 
the cos 20 asymmetry comes only from the Boer-Mulders 
effect in the region ^ Q^, and in this region the 
following relation hold: 



Xi 



X2 



(18) 



where y is the rapidity of the ^* /Z. We can also express 
the cross section of the Drell-Yan process depending on 
y and with an additional Jacobian determinant: 
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From Eq.Q, the azimuthal dependent terms are the 
second and the third terms with cos 2(j) and sin 2(j) forms 
respectively. However, the sin 2(f> term is suppressed, 



which can be found from ()T0|) and ([141) . shown in 
Ref. "4?!, we can write the coefficient of cos 20 term in 
Eq.(l5]) Waa into two parts, the perturbative QCD effect 
and the Boer-Mulders effect W^^. Then using 



pQCD 



AA 



an approximate Lam- Tung relation 2W^^^ — Wl ~ 0, 
one can give the cos 20 asymmetry caused by the Boer- 
Mulders effect: 
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Comparing (HJ and ([6]), and neglecting the Wl in the 
denominator because Wl <C Wt at low qx region, we 
can get the following relation: 
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In the numerical calculation, we choose the values of 
parameters in the Boer-Mulders function as those in 
Ref . l43l. 571. There is still an unsettled factor w which 



might be flavor dependent in the parametrization, be- 
cause it will be canceled in the product of two Boer- 
Mulders functions of quark and antiquark. It can range 
in the region 0.48 < w < 2.1, which is limited by the pos- 
itivity bounds 4^, 43, 58|. However, in the pp Drell-Yan 
process, it has the product of two Boer-Mulders functions 
of two quarks or two antiquarks which will not cancel the 
factor Lo. The cross section has different behavior with 
different values for oj. Therefore, we can learn additional 
information of the Boer-Mulders function from pp Drell- 
Yan processes. In this work, we choose three different 
values for a; = 0.5, lo = 1 and w = 2 to calculate 
and show their different behavior. 

In order to give with respect to a parameter y, Q 
or qx , we should integrate for the other parameters of the 
numerator and the denominator in Eg. (1211) respectively. 
The integral over need to be cut off &i qx = 2 GeV, 
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FIG. 1: The Q-dependent cos 20 azimuthal asymmetry caused 
by the Boer-Mulders effect in the unpolarized p+p — >■ 'y* /Z -\- 
X — >■ Z^/^ + X process at Z mass region. The solid, dashed, 
and dotted curves correspond to = 1, oj = 2, and oj = 0.5 
respectively. 
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FIG. 2: The rapidity y-dependent cos 20 azimuthal asym- 
metry caused by the Boer-Mulders effect in the unpolarized 
p+p — >■ 7*/Z-|-X — >■ +X process at Z mass region. The 
solid, dashed, and dotted curves correspond to w = 1, a; = 2, 
and u) — 0.5 respectively. 



because intrinsic transverse momentum plays a signifi- 
cant role at low qx and the fitting for the parameters 
has excluded the data with qr > 2 GeV. 



Comparing 
taken by CDF 



with the angular distribution form 



da 
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will contribute to ^2 caused by the Boer-Mulders 
effect at low qx- 

In summary, we calculated the cos20 azimuthal asym- 
metry in the unpolarized pp Drell-Yan dilepton produc- 
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FIG. 3: The gT-dependent cos2<^ azimuthal asymmetry 
caused by the Boer-Mulders effect in the unpolarized p + p — > 
7*/Z + X V^V -f X process at Z mass region. The solid, 
dashed, and dotted curves correspond to cj = 1, cj = 2, and 
= 0.5 respectively. 



tion processes in the Z mass region at CDF kinematic 
domain. It can be measured by experimental detection 
of the Lam- Tung relation violation. It is possible to study 
the spin structure of hadrons in unpolarized collision pro- 
cesses around Z mass region at Tevatron. In addition, 
the processes can give more significant information 
of the Boer-Mulders function than processes. It can 
help us to settle the factor to in the parametrization, and 
the prediction that the Boer-Mulders function have dif- 
ferent signs in SIDIS and Drell-Yan processes 21| also 
awaits experimental confirmation. Therefore the avail- 
able data of Tevatron are ideal to investigate the spin 
structure of nucleons via the unpolarized process at 
the Z pole. Besides, the GSTPANDA experiment [11] t^^, 
will run unpolarized Drell-Yan processes with colliding 
at s = 30 GeV^, and PAX experiment [60| may preform 
unpolarized Drell-Yan process with the fixed target 
mode at s = 45 GeV^. They will provide us an envi- 
ronment to study the Boer-Mulders effect at J/V" and 
T peaks and to understand the structure of nucleons. 
All these Drell-Yan experiments will give us signif- 
icant promotion in understanding the hadron structure 
and non-perturbative QCD properties. 
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